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Abstract 

A search for Bg oscillations is performed using a sample of semileptonic b-hadron 
decays collected by the ALEPH experiment during 1991-1995. Compared to previ- 
ous inclusive lepton analyses, the proper time resolution and b-flavour mistag rate 
are significantly improved. Additional sensitivity to Bg mixing is obtained by iden- 
tifying subsamples of events having a Bg purity which is higher than the average for 
the whole data sample. Unbinned maximum likelihood amplitude fits are performed 
to derive a lower limit of Am s > 9.5 ps -1 at 95% CL. Combining with the ALEPH 
D~ based analyses yields Am s > 9.6 ps _1 at 95% CL. 
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1 Introduction 



Flavour non-conservation in charged weak current interactions allows mixing between the 
B° and B° flavour states. The proper-time probability density function of a B° meson 
which is known to have mixed oscillates. The oscillation frequency is proportional to Am s , 
the mass difference between the mass eigenstates. Within the framework of the Standard 
Model, a measurement of the ratio Am s /Amd (Am^ being the mass difference in the 
B° — B[J system) would allow the extraction of the ratio of Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix elements |Vt s /Vtd|- 

Although the slower B° oscillations are now well established, the faster B° oscillations 
remain to be detected. Previous ALEPH analyses searching for B° oscillations have either 
been based on semi-exclusive selections in which a D~ is fully reconstructed [[[], |2| or on 
more inclusive lepton selections ||, f|, ||. Although the latter suffer from a lower B° purity 
and poorer proper time resolution they have the advantage of larger statistics. 

The analysis presented here is also based on an inclusive lepton sample. Compared to 
the previous ALEPH inclusive lepton analysis the following improvements are made 
to increase the sensitivity to B° mixing. 

• Decay length resolution: An improved decay length resolution is obtained by 
applying tight selection cuts to remove events likely to have misassigned tracks 
between the primary and the B° vertex. In addition an estimate of the decay 
length uncertainty is used on an event-by-event basis, rather than assuming the 
same average decay length uncertainty for all events, as used in previous analyses. 

• Boost resolution: A nucleated jet algorithm is used for an improved estimate of 
the momentum of the b-hadrons. 

• B° purity classes: Various properties of the events, such as the charge of the 
reconstructed b-hadron vertex and the presence of kaons are used to enhance the 
fraction of Bj? in subsamples of the data. 

• Initial and final state tagging: The b-flavour tagging method previously used for 
the D~ based analyses [l], 0] is applied. In this method discriminating variables are 
used to construct mistag probabilities and sample composition fractions estimated 
on an event-by-event basis. As a result, all events are tagged and the effective mistag 
rate is reduced. 



This paper details these improvements and is organized as follows. After a brief 
description of the ALEPH detector, the event selection is described in Section |^ and 
the B° purity classification procedure in Section |j. The next two sections explain the 
proper time reconstruction and the procedure for tagging the initial and final state b 
quark charge. The likelihood function is presented in Section |7| and the Am s results in 
Section ||. In Section | the systematic uncertainties are described, and in Section [10 



additional checks of the analysis presented. Finally the combination of this analysis with 



the ALEPH D s based analyses is described in Section |TT 
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2 The ALEPH detector 



The ALEPH detector and its performance from 1991 to 1995 are described in detail 
elsewhere || 0, and only a brief overview of the apparatus is given here. Surrounding 
the beam pipe, a high resolution vertex detector (VDET) consists of two layers of double- 
sided silicon microstrip detectors, positioned at average radii of 6.5 cm and 11.3 cm, 
and covering 85% and 69% of the solid angle respectively. The spatial resolution for the 
r(f) and z projections (transverse to and along the beam axis, respectively) is 12 /xm at 
normal incidence. The vertex detector is surrounded by a drift chamber with eight coaxial 
wire layers with an outer radius of 26 cm and by a time projection chamber (TPC) that 
measures up to 21 three-dimensional points per track at radii between 30 cm and 180 cm. 
These detectors are immersed in an axial magnetic field of 1.5 T and together measure 
the momenta of charged particles with a resolution u{p)/p = 6 x 10 -4 £>t © 0.005 (pr 
in GeV/c). The resolution of the three-dimensional impact parameter in the transverse 
and longitudinal view, for tracks having information from all tracking detectors and two 
VDET hits (a VDET "hit" being defined as having information from both r<p and z views), 
can be parametrized as o = 25 /jm + 95 fim/p {p in GeV/c). The TPC also provides up 
to 338 measurements of the specific ionization of a charged particle. In the following, 
the dE/dx information is considered as available if more than 50 samples are present. 
Particle identification is based on the dE/dx estimator Xw (xk), defined as the difference 
between the measured and expected ionization expressed in terms of standard deviations 
for the 7r (K) mass hypothesis. The TPC is surrounded by a lead/proportional-chamber 
electromagnetic calorimeter segmented into 0.9° x 0.9° projective towers and read out 
in three sections in depth, with energy resolution a(E)/E = 0.18/\AB + 0.009 (E in 
GeV). The iron return yoke of the magnet is instrumented with streamer tubes to form 
a hadron calorimeter, with a thickness of over 7 interaction lengths and is surrounded by 
two additional double-layers of streamer tubes to aid muon identification. An algorithm 
combines all these measurements to provide a determination of the energy flow |7j with 
an uncertainty on the total measurable energy of cr(E) = (0.6 J E /GeV + 0.6) GeV. 



3 Event selection 

This analysis uses approximately 4 million hadronic Z events recorded by the ALEPH 
detector from 1991 to 1995 at centre of mass energies close to the Z peak and selected with 
the charged particle requirements described in Ref. ||. It relies on Monte Carlo samples 
of fully simulated Z — > qq events. The Monte Carlo generator is based on JETSET 
7.4 || with updated branching ratios for heavy flavour decays. Monte Carlo events are 
reweighted to the physics parameters listed in Table [l]. 

Events for which the cosine of the angle between the thrust axis and the beam axis is 
less than 0.85 are selected. Using the plane perpendicular to the thrust axis, the event is 
split into two hemispheres. Electrons and muons are identified using the standard ALEPH 
lepton selection criteria ||10|| . Events containing at least one such lepton with momentum 
above 3 GeV/c are kept. The leptons are then associated to their closest jet (constructed 
using the JADE algorithm |TT| with y cnt = 0.004) and a transverse momentum pt with 



respect to the jet is calculated with the lepton momentum removed from the jet. Only 
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leptons with pr > 1.25 GeV/c are selected. In the case that more than one lepton in an 
event satisfies this requirement, only the lepton with the highest momentum is used as a 
candidate for a B° decay product. 

The e + e~ interaction point is reconstructed on an event-by-event basis using the 
constraint of the average beam spot position and envelope JT2] . 



A charm vertex is then reconstructed in the lepton hemisphere using the algorithm 
described in Ref . || . Charged particles in this hemisphere (other than the selected lepton) 
are assigned to either the interaction point or a single displaced secondary vertex. A 
three-dimensional grid search is performed for the secondary vertex position to find the 
combination of assignments that has the greatest reduction in x 2 as compared to the case 
when all tracks are assumed to come from the interaction point. Tracks are required to 
come within 3a of their assigned vertex. The position resolution of this "charm vertex" 
is subsequently improved by removing those tracks having a momentum below 1.5 GeV/c 
or an impact parameter significance relative to the charm vertex larger than 1. 4a. The 
remaining tracks are then re-vertexed to form the reconstructed "charm particle" . If only 
one track passes the requirements, it serves as the charm particle. The event is rejected 
if no track remains, or none of the tracks assigned to the charm vertex have at least 
one vertex detector hit. The charm particle is then combined with the lepton to form a 
candidate b-hadron vertex. The lepton is required to have at least one vertex detector 
hit and the x 2 P er degree of freedom of the reconstructed b-hadron vertex is required to 
be less than 25. 

The energy E c of the charm particle is estimated by clustering a jet, using the JADE 
algorithm, around the charged tracks at the charm vertex until a mass of 2.7 GeV/c 2 is 
reached. To reduce the influence of fragmentation particles on the estimate of E c , charged 
and neutral particles with energies less than 0.5 GeV are excluded from the clustering |l6fl . 
The neutrino energy E u is estimated from the missing energy in the lepton hemisphere 



taking into account the measured mass in each hemisphere |i7fl . Assuming the direction 
of flight of the b-hadron to be that of its associated jet, an estimate of the b-hadron 
mass can be calculated from the energy of the neutrino and the four-vectors of the charm 
particle and the lepton. 

In order to improve the rejection of non-b background or b events with a badly 
estimated decay length error, the following additional cuts are applied 



• the momentum of the charm particle must be larger than 4 GeV/c; this cut is 
increased to 8 GeV/c when the angle between the charm particle and the lepton is 
less than 10°; 

• the reconstructed mass of the b-hadron must be less than 8 GeV/c 2 ; 

• the missing energy in the lepton hemisphere must be larger than —2 GeV; 

• the angle between the charm particle and the lepton must be between 5° and 30°; 

• the angle between the charm particle and the jet must be less than 20°. 

Although the total efficiency of these additional requirements is 35%, the average decay 
length resolution of the remaining events is improved by a factor of 2 and the amount of 
non-b background in the sample reduced by a factor close to 4. In addition the average 



3 



Table 1: Values of the physics parameters assumed in this analysis. 



Physics parameter 


Value and uncertainty 


Reference 


B° lifetime 
B° lifetime 
B + lifetime 
b-baryon lifetime 
Amd 


1.49 ±0.06 ps 
1.57 ±0.04 ps 
1.67 ±0.04 ps 
1.22 ±0.05 ps 
0.463 ±0.018 ps" 1 


[13 

[13] 

[13] 

(til 
@ 


R h = B(Z -> bb)/B(Z -> qq) 
i? c = B(Z -> cc)/B(Z -> qq) 

/bo = B(b -> B s °) 
/bo =/ B+ = i3(b^B°,B + ) 
/b-baryon = B(b -> b-baryon) 


0.2170 ±0.0009 
0.1733 ±0.0048 

_ „ _ _i_n ni k 

u.iuo_ 015 
u.dyo_ 002 o 
0.1061^1? 


B 

[0 
[0 



B(b -> t) 
B(b -> c -> £) 
B(b -> c -> f) 

B(c -> £) 


0.1112 ±0.0020 
0.0803 ±0.0033 
0.0013 ±0.0005 
0.098 ±0.005 










0.702 ±0.008 






Table 2: Lepton candidate sources (%), as estimated from Monte Carlo. Quoted uncertainties 
are statistical only. 





R° 


other b-hadrons 


charm 


uds 


10.35 ±0.08 


38.53 ±0.13 


47.86 ±0.14 


2.31 ±0.06 


0.95 ±0.05 



momentum resolution of the sample is significantly improved. A total of 33023 events 
survive after all cuts. 

4 Bg purity classes 

Table [| shows the composition of the final event sample obtained assuming the physics 
parameters listed in Table [I] and reconstruction efficiencies determined from Monte Carlo. 
The average B° purity in the sample is estimated to be 10.35%. 

The sensitivity of the analysis to B° mixing is increased by splitting the data into 
subsamples with a B° purity larger or smaller than the average and then making use 
of this knowledge in the likelihood fit. Classes are constructed based on (i) the track 
multiplicity at the charm vertex, (ii) the number of identified kaon candidates and (iii) the 
charge correlation between the tracks at the charm vertex and the lepton. The definition 
of the eleven classes used in this analysis is given in Table [| As the last class contains 
those events which do not satisfy the criteria of the preceding classes, the classification 
procedure does not reject events. 

For an odd (even) number of charged tracks assigned to the charm vertex, the 
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Table 3: Definition of the eleven B[? purity classes. Column 1 gives the number of charged 
tracks at the charm vertex. Column 2 shows whether the charge of these tracks are the same (S) 
or opposite (O) to that of the lepton, the tracks being ranked in order of decreasing momentum. 
Column 3 indicates the subclasses based on the presence of kaon or cj> candidates at the charm 
vertex. Column 4 shows the fraction of data events in each class. Column 5 gives the Bg purity 
in each class, as estimated from Monte Carlo. Quoted uncertainties are statistical only. 



Number 
of tracks 


Charge 
correlation 


Kaon 
requirements 


Traction m 

data (%) 


B° purity (%) 


1 





1 kaon 
kaon 


3.8 
14.9 


24.0 ± 0.6 
14.7 ± 0.3 


2 


OS, SO 
OS, SO 
OS, SO 
OS, SO 




kaon 

1 kaon 

2 kaons 


1.2 
17.8 
17.4 

2.3 


21.1 ± 1.0 
7.0 ± 0.2 
5.2 ± 0.1 
8.4 ± 0.5 


00 




8.3 


16.7 ± 0.4 


3 


OOS 
OSO 
SOO 




2.9 
3.8 
3.9 


19.4 ± 0.6 
18.0 ± 0.5 

14.5 ± 0.5 


remainder 


23.6 


5.7 ± 0.1 



reconstructed charge of the b-hadron vertex is more likely to be zero (non-zero), and 
therefore the probability for the hemisphere to contain a neutral b-hadron is enhanced 
(reduced). For events having two oppositely charged tracks at the charm vertex, the Bg 
purity is 6.7%, which is lower than the average purity. For this large subsample of events, 
the presence of kaon candidates and consistency with the mass are used to recover some 
sensitivity to the Bg. In this procedure, kaon candidates are defined as charged tracks 
with momentum above 2 GeV/c satisfying Xn + Xk < and \xk\ < 2, and a candidate 
is defined as a pair of oppositely charged tracks with an invariant mass between 1.01 and 
1.03 GeV/c 2 (assuming kaon masses for the two tracks). 

Monte Carlo studies indicate that this classification procedure is effectively equivalent 
to increasing the statistics of the sample by 28%. 

5 Proper time reconstruction and resolution 

An estimate, I, of the decay length of each b-hadron candidate is calculated as the distance 
from the interaction point to the b-hadron vertex projected onto the direction of the 
jet associated to the lepton. This decay length includes a global correction of —78 /im, 
determined using Monte Carlo events. This small offset is due to the vertex reconstruction 
algorithm, which assumes that all lepton candidates in b events come from direct b —>■ £ 
decays. Figure [TJa shows the Monte Carlo distribution oil — lo for b events, where Iq is the 
true decay length. An event-by-event decay length uncertainty, 07, is estimated from the 
covariance matrices of the tracks attached to the vertices. This can be compared with the 
true error, (I — l ), by constructing the pull distribution, (I — Zo)/°7- A fit to this Monte 
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1-1 ( cm > (g-8oVgo 

Figure 1: Decay length resolution (a) and relative boost term resolution (b) for all b-hadrons; 
the curves are the result of fits of the sum of two Gaussian functions with relative fractions and 
widths as indicated. 



Table 4: Double-Gaussian parametrizations 
resolution obtained from Monte Carlo. 



Parametrization of (I — Zq) /cz 


a 


Fraction /" 


Sigma Sf 


1 


0.849 ± 0.003 


1.333 ±0.005 


2 


0.151 ±0.003 


4.365 ± 0.033 



the decay length pull and relative boost term 



p 


3 arametrization 
Fraction ffj 


of (g - go) /g 
Sigma S@ 


1 

2 


0.723 ± 0.004 
0.277 ±0.004 


0.0713 ±0.0003 
0.2101 ±0.0012 



Carlo distribution of the sum of two Gaussian functions (a = 1, 2) yields the fractions, 
and sigmas, Sf , indicated in Table |j. These parameters are used to describe the observed 
tails when constructing the resolution function. 

The true boost term is denned as go = t /l , where t is the true proper time. An 
estimate of the boost term is formed using g = mB/pB±0.36 ps/cm. The average b-hadron 
mass, 771b, is assumed to be 5.3 GeV/c 2 and the reconstructed momentum is calculated 
as pb = J (E c + E v + Ei) 2 — m| where Ei is the measured lepton energy. The constant 
term is an average offset correction determined using Monte Carlo events; this results from 
the choice of the mass cut-off used in the nucleated jet algorithm described in Section [3|, 
which optimizes the relative boost term resolution. The distribution of (g — go)/ go, shown 
in Fig. |I|b, is parametrized with the sum of two Gaussian functions; Table |] shows the 
corresponding fractions, f@, and sigmas, S^, determined with Monte Carlo events. 

The proper time of each b-hadron candidate is computed from the estimated decay 
length and boost term as 

t = lg, (1) 
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Figure 2: The proper time resolution for b events in various intervals of true proper time to 
(in ps). The curves display the corresponding resolution assumed in the likelihood as obtained 
from Eq. The RMS values indicated, are derived from the data points shown. 



and its proper time resolution function is parametrized with the sum of four Gaussian 
components, 



2 2 



Res(Mo) = yy f?rt-i= 



cxp 



1 t-U 



2 U^(t ) 



(2) 



where ff = ff^ff and f}' = 1 — ff, and where the event-by-event resolution a a(3 of each 
component, given by 



*°"(*o) = MgSFSfa) +(toS**Si 



S), 



(3) 
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Proper time (ps) 



Figure 3: The reconstructed proper time distributions of the selected events in data. The 
contributions from the various components are indicated. The curve is the result of the fit 
described in Section |5[ 



includes the explicit dependence on t . This parametrization implicitly assumes that any 
correlation between the decay length resolution and the relative boost resolution is small, 
as confirmed by Monte Carlo studies. 

The scale factors Sf at and / dat are introduced to account for a possible discrepancy 
between data and Monte Carlo, both in the amount of tail in the decay length pull 
(// dat ) an d in the estimate of <j\ itself (S dat ). In a similar fashion, the inclusion of the 
parameter Sg 3 * 2 allows possible systematic uncertainties due to the boost resolution to be 
studied. By definition all these factors are set to unity when describing the resolution of 
simulated events. Figure |2| shows, for various intervals of true proper time, the proper 
time resolution of simulated b events together with the parametrization obtained from 
Eq. (fj). The parametrization is satisfactory, especially for small proper times. 

In order to measure Sf at and / dat in the data, a fit is performed to the reconstructed 
proper time distribution of the selected sample of real events. This is performed using the 
likelihood function described in Section [7], modified to ignore tagging information. Fixing 
all physics parameters to their central values given in Table [j], the likelihood is maximized 
with respect to Sf at and / dat . The fit reproduces well the negative tail of the proper time 
distribution (see Fig. |3p, showing that the resolution is satisfactorily described by the 
chosen parametrization. The fitted values Sf at = 1.02±0.03 and / dat = 1.20±0. 09 indicate 
that the decay length resolution in the data is somewhat worse than that suggested by 
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Figure 4: Schematic drawing indicating the initial and final state tags used in this analysis. 



the Monte Carlo simulation. 



6 Initial and final state tagging 

The flavour state of the decaying B° candidate is estimated from the charge of the 
reconstructed lepton. This final state tag is incorrect if the lepton is from the b — > c — > £ 
decay (6.1% of the b events in the sample) as the charge of the lepton is reversed. The 
flavour state at production is estimated using three initial state tags. A B° candidate 
is "tagged as unmixed (mixed)" when the reconstructed initial and final flavour states 
are the same (different). By definition, candidates from charm, uds, or non-oscillating 
b-hadron backgrounds are correctly tagged if they are tagged as unmixed. 

The tagging power is enhanced by the means of discriminating variables which have 
some ability to distinguish between correctly tagged and incorrectly tagged candidates. 
This approach was first used in the ALEPH D~-lepton analysis and refined for the 
DjT-hadron analysis ||. In contrast to what was performed in Refs. |T| and 0, an event 
is considered to be mistagged if either the initial or final state is incorrectly tagged, but 
not both. 

For each B° candidate, one of the tags described below is used to determine the initial 
state (see also Fig. |j). 

• Opposite lepton tag: Leptons with momentum larger than 3 GeV/c are searched 
for in the hemisphere opposite to the B° candidate. The sign of the lepton with the 
highest transverse momentum pr(^o) tags the nature of the initial b quark in the 
opposite hemisphere. It takes precedence over the other tags if it is available. 

• Fragmentation kaon tag: The fragmentation kaon candidate is defined as the 
highest momentum charged track within 45° of the B° direction, identified, using 



9 



the vertexing algorithm described in Section 2, as being more likely to come from the 
interaction point than the charm vertex, and satisfying xk < 0.5 and xk — Xw > 0.5. 

The sign of the fragmentation kaon candidate tags the sign of the b quark in the 
same hemisphere. It is used if no opposite hemisphere lepton tag is found. 

• Opposite hemisphere charge tag: The opposite hemisphere charge is defined as 

oppo 

o oppo ) (4) 

E Mr 

i 

where the sum is over all charged particles in the opposite hemisphere, pi is the 
momentum of the i th track projected on the thrust axis, its charge and k = 0.5. 
The sign of Q D tags the initial state of the b quark in the opposite hemisphere. This 
tag is always available but has the largest mistag probability of the three tags. It is 
used only if no other tag is available. 

The events are sorted into five exclusive classes based on the availability and results 
of the three tags. The definition of these tagging classes and the list of the discriminating 
variables associated with each class are given in Table 0. The variable Q s is the sum 
of the charges of all the tracks in the same hemisphere and carries information on the 
initial state of the B°. As the sum of charges of tracks originating from the decay of a 
neutral particle is zero, it is independent of whether the B° decays as a B° or a B°. The 
variable Z K is defined as Z K = p K / (E beam — E B ), where pk is the kaon momentum, E heiLm 
the beam energy and the B° candidate energy. The inclusion of the reconstructed B° 
proper time t takes into account that the mistag probability of the fragmentation kaon 
tag increases as the B° vertex approaches the primary vertex, due to the misassignment of 
tracks between the primary and secondary vertices. The use, for all classes, of the variable 
Pt{£ s ), the transverse momentum of the lepton from the B° candidate decay, reduces the 
deleterious effect of b — > c — ► I on the final state mistag. 

The mistag probability, rj, for the B° signal events in each class, as well as the 
probability distributions of each discriminating variable Xi for correctly and incorrectly 
tagged signal events, r^Xi) and Wi(xi), are estimated from Monte Carlo. The various 
discriminating variables chosen in each class, x±,X2, ■ ■ ., are combined into a single effective 
discriminating variable x eS , according to the prescription developed for the D~ based 
analyses p], |2j . This new variable is defined as 

x eS = VW 1 (X 1 )W 2 (X 2 ) ■■■ _ 

{l-r])ri{xi)r 2 {x 2 ) ■■■ +r]Wi{x 1 )w 2 {x 2 ) 

and takes values between and 1. A small value indicates that the B° oscillation is likely 
to have been correctly tagged. 

To allow use of the discriminating variables in the likelihood fit, the probability density 
functions Gj kl (x eS ) of x eS are determined for each lepton source j, in each tagging class k 
and in each B° purity class I, separately for the correctly (c = +1) and incorrectly (c = —1) 
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Table 5: The tag and discriminating variables used in each class. The quantities S(Q Q ), S(K) 
and S(£ ) are the signs of the opposite hemisphere charge, the fragmentation kaon and the 
opposite side lepton. Classes 3-5 all use the sign of the opposite hemisphere lepton as the initial 
state tag. For Class 3 no fragmentation kaon candidate is identified. For Class 4 (Class 5) a 
fragmentation kaon candidate is found whose charge is the same as (opposite to) the charge 
of the opposite hemisphere lepton. Purity and mistag rates are estimated from Monte Carlo. 
Quoted uncertainties are statistical only. 
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pt(£ ) 
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pt(£ ) 
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t 




Pt{Q 


Pt(Q 


Pt(Q 


Pt(Q 


pt(£ s ) 


Fraction in data (%) 


71.4 ± 0.2 


11.9 ± 0.2 


14.2 ± 0.2 


1.3 ± 0.1 


1.2 ± 0.1 


B s ° purity (%) 


9.8 ± 0.1 


13.1 ± 0.3 


10.1 ± 0.2 


15.6 ± 1.0 


11.8 ± 0.8 


B° mistag (%) 


38.6 ± 0.5 


28.9 ± 1.0 


34.0 ± 1.1 


16.1 ± 2.3 


55.9 ± 3.5 


B° effective mistag (%) 


32.4 


24.0 


24.5 


12.5 


22.3 


B° d mistag (%) 


38.4 ± 0.2 


48.5 ± 0.7 


35.4 ± 0.5 


35.5 ± 2.0 


39.9 ± 2.0 


other B mistag (%) 


37.6 ± 0.2 


61.4 ± 0.5 


34.2 ± 0.5 


43.8 ± 1.8 


24.1 ± 1.4 


charm mistag (%) 


38.2 ± 1.4 


54.2 ± 3.2 


14.2 ± 3.1 


50.0 ± 50.0 


8.6 ± 8.2 


uds mistag (%) 


47.8 ± 2.8 


56.9 ± 6.0 


46.0 ± 12.9 


50.0 ± 50.0 


50.0 ± 50.0 



tagged events. This determination (as well as the estimation of the corresponding mistag 
probabilities rjjki) is based on Monte Carlo. 

The enhancement of the tagging power provided by the variable, x eS , depends on the 
difference between the G^ kl {x cS ) and G~ kl (x eS ) distributions, and can be quantified in 
terms of effective mistag rates, as described in Ref. ffl. The effective mistag rates for the 
B° signal in the five tagging classes are given in Table [|. This table also indicates B° 
purity and the mistags for all background components. The overall average B° effective 
mistag is 29%. 

Figure |5] displays the distribution of x eS in each of the tagging classes; a good agreement 
is observed between data and Monte Carlo. The systematic uncertainties associated with 
the tagging procedure are considered in Section |. 
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Figure 5: Distribution of x in each tagging class in data (points) and Monte Carlo (histogram). 



7 Likelihood function 

Each b-hadron source has a different probability distribution function for the true proper 
time, to, and for the discrete variable, /io, defined to take the value —1 for the mixed case 
or +1 for the unmixed case. Assuming CP conservation and equal decay widths for the 
two CP eigenstates in each neutral B-meson system, the joint probability distribution of 
t and [io can be written as 

e ~to/Tj 

PjiHo, t ) = — [1 + /io cos (Arrij t )} , (6) 

where Tj and Arrij are the lifetime and oscillation frequency of b-hadron source j (with 
the convention that Arrij = for non-oscillating b-hadrons). 

The efficiency for reconstructing the b-hadron vertex depends on the true proper time. 
The stringent selection cuts described in Section |3| are designed to reduce the fraction of 
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fragmentation tracks assigned to the charm vertex, consequently causing a loss of efficiency 
at small proper times. Similarly at large proper times the efficiency also decreases as one 
is less likely to include a fragmentation track at the charm vertex and therefore more 
likely to fail the requirement of the charm vertex being assigned at least one track. The 
efficiencies €j{to) are parametrized separately for each b-hadron component j. They are 
independent of whether the b-hadron candidate is tagged as mixed or unmixed. 

The joint probability distribution of the reconstructed proper time t and of /io is 
obtained as the convolution of Pj((J>o,t ) with the event-by-event resolution function 
Res(t,t ) (Section |5|) and takes into account the observed dependence of the selection 
efficiency on true proper time: 



ej (t )pj (/i , t )Res(t, t ) dt 

hj(tM>,t) = ~ 755 1 • (7) 

' e i (to)-e-* o/T ^t 



For the lighter quark backgrounds, hj(—l,t) = as these sources are unmixed 
by definition, and hj(+l,t) are the reconstructed proper time distributions. These 
distributions are determined from Monte Carlo samples and are parametrized as the sum 
of three Gaussian functions. 

The likelihood function used in this analysis is based on the values taken by three 
different variables in the selected data events. These variables are the reconstructed 
proper time t, the tagging result /x, taking the value —1 for events tagged as mixed or 
+ 1 for those tagged as unmixed, and the effective discriminating variable x eS . The use 
of the discriminating variable x eS in the likelihood function is reduced to the use of two 
sets of functions of x eS , Xj k i(x eS ) and Yjki(x eS ) (described below), whose values can be 
interpreted as event-by-event mistag probabilities and fractions of the different lepton 
sources respectively. The likelihood of the total sample is written as 

11 purity 5 tagging N^i events 

£=c n n n /«(*£. z^,**), w 

l k i 

where C is a constant independent of b-hadron oscillation frequencies and lifetimes, N k i 
is the number of selected candidates from B° purity class / falling in tagging class k, and 
where 

5 sources 

f kl (x eS ,^t) = J2 ^H^ cff )[(l-X jfc Kx eff ))/ ij (/i,t)+X^K^ fr )^(-/i,t)] (9) 
j 

sums over the 5 different lepton sources considered to comprise the sample (see Table 0) . 
The event-by-event quantities Xjki(x ) and Yj k i(x ) are computed from the distributions 
Gj kl (x eS ) and mistag probabilities r]j k[ introduced in Section ^ 

X jki (x ) - rj jM , Y jk i(x )-a jk i—— — - — -— ^ , (1UJ 

where Gj k i(x cS ) = (l—rij k i)G~j" kl ( y x eS )+rij k iGj kl (x eS ) and where ctjki are the source fractions, 

.5 

j-- 



satisfying ^1^^1 = 1- 
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Figure 6: Negative log-likelihood difference with respect to the minimum as a function of Am s . 



8 Results for Am s 

Assuming the values for the physics parameters given in Table fj, the variation in the data 
of the log-likelihood, as a function of the free parameter Aro s , is shown in Fig. [|. The 
difference in log-likelihood is plotted relative to its global minimum and remains constant 
for Am s larger than 20 ps -1 . The global minimum occurs at Am s = 15.9 ± 1.6(stat.) ps -1 
but is not sufficiently deep to claim a measurement. 

In order to extract a lower limit on Am s and to facilitate combination with other 
analyses, the results are also presented in terms of the "amplitude" fit. In this method jnj 
the magnitude of B° oscillations is measured at fixed values of the frequency Am s , using 
a modified likelihood function that depends on a new parameter, the B° oscillation 
amplitude A. This is achieved by replacing the probability density function of the B° 
source given in Eq. (H) with 

g-*o/T s 

— [1 + (iqA cos (Am s t )] . (11) 

2r s 

For each value of Am s , the new negative log-likelihood is then minimized with respect to 
A, leaving all other parameters (including Am s ) fixed. The minimum is well behaved and 
very close to parabolic. At each value of Am s one can thus obtain a measurement of the 
amplitude with Gaussian error, A ± cr^ at . If Am s is close to the true value, one expects 
A = 1 within the estimated uncertainty; however, if Am s is far from its true value, a 
measurement consistent with A = is expected. 

The amplitude fit results are displayed in Fig. [?] as a function of Am s . A peak in the 
amplitude, corresponding to the minimum observed in the negative log-likelihood, can be 
seen around Am s = 16 ps -1 . At this value, the measured amplitude is 2.2 a away from 
zero; as for the likelihood, this is not significant enough to claim a measurement of Am s . 
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Figure 7: Measured Bg oscillation amplitude as a function of Am s for this analysis. The error 
bars represent the la total uncertainties, and the shaded bands show the one-sided 95% CL 
contour, with and without systematic effects included. 



A value of Am s can be excluded at 95% CL if ^4+ 1.645 04 < 1. Taking into account 
all systematic uncertainties described in the next section, all values of Am s below 9.5 ps _1 
are excluded at 95% CL. The sensitivity, estimated from the data as the value of Am s 
at which 1.645 04 = 1, is 9.6 ps -1 . Ignoring systematic uncertainties would increase the 
95% CL lower limit and sensitivity by 0.1 ps -1 and 0.6 ps -1 respectively. 
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9 Systematic uncertainties 



The systematic uncertainties on the B° oscillation amplitude a s J st are calculated, using 
the prescription of Ref . JL9| , as 



a syst = ^new _ ^nom + ( X _ J"m)^_^_ ^ 

where the superscript "nom" refers to the amplitude values and statistical uncertainties 
obtained using the nominal values for the various parameters and "new" refers to the new 
amplitude values obtained when a single parameter is changed and the analysis repeated 
(including a re-evaluation of the distributions of the discriminating variables used for the 
b-flavour tagging). The total systematic uncertainty is the quadratic sum of the following 
contributions. 

• Sample composition: The systematic uncertainty on the sample composition is 
obtained by varying the assumed values for the b-hadron fractions / B o, /b-baryon and 
the various lepton sources (b — > £, b — > c — > £, etc . . . ) by the uncertainties quoted 
in Table [TJ. The statistical uncertainty on the purities determined from Monte Carlo 
is also propagated. 

A comparison of data and Monte Carlo fractions for the different B° purity classes 
shows small deviations, the largest relative difference of 16% occurring in the first 
class of Table §. The systematic uncertainty due to the B° purity classification 
procedure is evaluated by shifting, in each class, all five purities (B° , B°, . . . ) in 
the direction of their respective overall averages, Hj given in Table [|, by a fraction 
7 = ±20% of their differences with respect to these averages: 

atjki — > ocjki + l(a jM - Uj) . (13) 

As this is performed coherently in all B° purity classes, the procedure is rather 
conservative and ensures that the overall average purities remain unchanged. Not 
using the B° purity classification would decrease the Aro s statistical sensitivity by 
0.7 ps- 1 . 

For the fraction of charm and uds backgrounds a relative variation of ±25% 
is propagated, as suggested from a comparison between data and Monte Carlo 
performed in Ref. [2(| . 



Proper time resolution: For the systematic uncertainty on the proper time 
resolution, the correction factors presented in Tables 4 and 5 are varied by ±ler. 
The scale factors (S, dat = 1.02 ± 0.03 and / ; dat = 1.20 ± 0.09) for the decay length 
resolution, obtained from the lifetime fit to the data, are also varied by their 
measured uncertainty. In addition, a possible bias of ±0.055 ps/cm is considered 
on the determination of the boost term; this value corresponds to the observed 
shift between the measured and simulated boost term distributions and represents 
approximately 1% of the average boost term. Finally the boost term resolution is 



the close agreement between the measured and simulated boost distributions. 



given a relative variation of ±10% (Sg = 1.0 ± 0.1), which is conservative given 
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• b-quark fragmentation: The average fraction of energy taken by a b-hadron 
during the fragmentation process, (Xe) = 0.702 ± 0.008, is varied by its measured 
uncertainty. The corresponding effects on the sample composition, mistags and 
resolutions are propagated. 

• Mistag: Based on data/Monte Carlo comparisons of the tagging variables, 
performed for the D~-based analyses [JJ, 0, absolute variations of ±0.8% for the 
first tagging class (opposite hemisphere charge) and ±2% for all other classes 
(fragmentation kaon and opposite lepton) are applied to the mistag rates. In 
addition, the ±ler statistical uncertainty from Monte Carlo is propagated. 

The changes in mistag due to variation of the b — > c — > £ fraction are included as 
part of the sample composition systematic uncertainty. 

• Lifetimes, Am,j, -Rb and R c : The values assumed for the various b-hadron 
lifetimes, Am<j, -Rb and R c are varied within the uncertainties quoted in Table [I]. 

• Difference in decay width: A possible decay width difference Ar s /T s between 
the two mass eigenstates of the B° meson has been ignored in the likelihood fit. 
The fit is therefore repeated with a modified likelihood assuming Ar s /T s = 0.27, 
equal to the theoretical prediction of Ref. |2"T|] , Ar s /r s = 0.16to]o9) P ms its quoted 
positive uncertainty. 

• Cascade bias: In the likelihood expression of Eq. each b-hadron component is 
treated using a single resolution function and mistag. No attempt is made to treat 
separately the b — > £ (direct) and b — * c — > £ (cascade) decays. While the former is 
characterized by a good proper time resolution and mistag, the latter has a degraded 
decay length resolution and a somewhat biased decay length because of the charm 
lifetime. In addition, the sign of the lepton is changed, leading to a different total 
mistag. To study the possible bias arising from the correlation between the poor 
decay length resolution and degraded tagging performance of the cascade events, 
two different fast Monte Carlo experiments are generated with a true value of 
Am s equal to 50 ps _1 . In the first the b-hadron decays are generated using the 
average mistag and resolution; in the second, the primary and cascade components 
are generated separately, each with their appropriate mistag and resolution. For 
both experiments, the corresponding amplitude plot is obtained using the likelihood 
described in Section ^, i.e. with average mistags and resolutions. 

The fast Monte Carlo experiment generated using the average b-hadron properties, 
yields an amplitude spectrum consistent with zero, as expected (since the fitting 
function is based on the same probability distributions as the fast Monte Carlo 
generator). In contrast, the experiment in which the direct and cascade decays are 
generated separately shows a small amplitude bias at low and very large Am s . Since 
the bias is small, especially in the region where the limit is set, and would cause the 
limit and sensitivity to be slightly underestimated, no attempt is made to correct for 
this effect; instead the deviations of the amplitude from zero observed are treated 
as a systematic uncertainty. 

The relative importance of the various systematic uncertainties, as a function of Am s , 
is shown in Table M. Except at low Am s the systematic uncertainties are generally 
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Table 6: Measurement of the B° oscillation amplitude, A, for various oscillation frequencies 
together with the statistical uncertainty, <r^ at , and the total systematic uncertainty, cr s J st ; a 
breakdown of a s J st in several categories of systematic effects is also given. 



Am s 


ps- 1 


5 ps" 1 


10 ps" 1 


15 ps" 1 


A 


-0.030 


-0.065 


0.303 


2.291 


_stat 
°A 


±0.099 


±0.267 


±0.590 


±1.271 


syst 

°A 


+0.340 
-0.340 


+0.223 
-0.235 


+0.232 
-0.324 


+0.801 
-0.582 


Systematic contributions: 










t~) r> 

- -Kb, -Kc 


+0.001 
-0.001 


+0.002 
-0.001 


+0.001 
-0.002 


+0.001 
-0.005 


- / B o = B{b -> BJ 


+0.046 
-0.035 


+0.146 
-0.112 


+0.133 
-0.109 


+0.217 
-0.173 


- /b-baryon = #(b -> b-baryon) 


+0.008 
-0.010 


+0.026 
-0.018 


+0.028 
-0.023 


+0.007 
-0.002 


- cnarm traction 


+0.012 
-0.012 


+0.019 
-0.016 


+0.021 
-0.018 


+0.051 
-0.043 


- uds traction 


+0.008 
-0.008 


+0.023 
-0.026 


+0.032 
-0.038 


+0.078 
-0.091 


- b — > -e, b— > c — > -e, b— > c — > -e, c — ► £ 


+0.065 
-0.013 


+0.000 
-0.055 


+0.000 
-0.121 


+0.464 
-0.000 


- purities (MC stat.) 


+0.047 
-0.041 


+0.078 
-0.070 


+0.076 
-0.075 


+0.104 
-0.108 


- B° purity classes 


+0.017 
-0.009 


+0.000 
-0.007 


+0.010 
-0.018 


+0.140 
-0.187 




+0.037 
-0.037 


+0.002 
-0.002 


+0.001 
-0.001 


+0.000 
-0.003 


- b-hadron lifetimes 


+0.033 
-0.000 


+0.000 
-0.046 


+0.027 
-0.037 


+0.282 
-0.000 


- decay length resolution 


+0.000 
-0.000 


+0.025 
-0.025 


+0.054 
-0.057 


+0.050 
-0.021 


- boost term resolution 


+0.010 
-0.010 


+0.030 
-0.033 


+0.048 
-0.059 


+0.205 
-0.191 


- b-fragmentation 


+0.023 
-0.000 


+0.012 
-0.070 


+0.067 
-0.085 


+0.509 
-0.403 


- b-flavour tagging 


+0.317 
-0.332 


+0.138 
-0.132 


+0.132 
-0.207 


+0.233 
-0.219 


- Aiyr s 


+0.000 
-0.002 


+0.012 
-0.000 


+0.011 
-0.000 


+0.018 
-0.000 


- cascade bias 


+0.060 
-0.000 


+0.000 
-0.087 


+0.000 
-0.085 


+0.000 
-0.069 



small compared to the statistical uncertainty. At Am s = 10 ps 1 , the most important 
contributions are from / B o and the b-flavour tagging. 

10 Checks 

Using a fast Monte Carlo generator which takes into account all details of the sample 
composition, the resolution functions, the mistag rates and the distributions of x cS , the 
average amplitude over many fast Monte Carlo experiments is found to be consistent 
with unity for Am s = Am^™ and with zero for any value of Am s if Am lrue = oo. The 
estimate, a^ at , of the statistical uncertainty on the amplitude has also been verified by 
studying the distribution of A/a^ for cases where A = is expected. The mean value 
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Figure 8: Measured Bj? oscillation amplitude as a function of Am s in the Z — > qq Monte Carlo. 
The error bars represent the la statistical uncertainties, the solid curve the one-sided 95% CL 
contour (systematic effects not included). The dotted line is 1.645a. The generated value of 
Am s is 3.33 ps -1 . 



and RMS of such a distribution obtained with fast Monte Carlo experiments generated 
with Am* rue = oo are found to be consistent with and 1. 

A likelihood fit for Am s performed on a Z — > qq Monte Carlo sample having the 
same statistics as the data and generated with a true value of Am s of 3.33 ps -1 yields 
Am s = 3.31 ± 0.12(stat.) ps -1 , in agreement with the input value. Performing an 
amplitude fit on the same Monte Carlo events yields the results shown in Fig. [8]. As 
expected, the amplitude is consistent with 1 at the true value of Am s . The sensitivity 
estimated from this Monte Carlo sample (ignoring systematic uncertainties) is 10.6 ps -1 , 
a little higher than that obtained from the data, 10.2 ps -1 , due to the slightly better 
decay length resolution in Monte Carlo. 

As a further check of the assumed mistags and sample composition, the analysis is 
used to measure Amd in the data. Fixing Aro s to 50 ps -1 and minimizing the negative 
log-likelihood with respect to Amd gives Amd = 0.451 ± 0.024(stat.) ps -1 , consistent 



with the latest world average of 0.463 ± 0.018 ps [fL3|1 . Figure |9] shows that the fitted 
B° oscillation amplitude is consistent with that observed in the Z — > qq Monte Carlo and 
has the expected value of 1 at the minimum of the negative log-likelihood. To check that 
the sample composition and mistags assumed for each B° purity class and tagging class 
are reasonable, a fit for the B° oscillation amplitude is performed separately in each class. 
At Amd = 0.451 ps~ x a value of A consistent with 1 is found in all classes; the largest 
deviation being 1.5cr st at in the last B° purity class ("remainder"). 
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Am d (ps ) 




Am d (ps ) 



Figure 9: Measured B° oscillation amplitude as a function of Am^ in (a) the data and (b) the 
2 — ► qq Monte Carlo. The error bars represent the la total uncertainties and the curves the 
one-sided 95% CL contour (systematic effects not included). 



11 Combination with D s analyses 

The amplitudes measured in this analysis and in the two ALEPH D~ analyses ffl, Q are 
combined. The small number of events common to both this analysis and the D~-lepton 
analysis are removed from the inclusive lepton sample before combining the results. The 
following sources of systematic uncertainty are treated as fully correlated: the values 
assumed for / B o, /b-baryon, Amd and the various b-hadron lifetimes, the b fragmentation, 
the decay length resolution bias in the Monte Carlo simulation Sf &t and / ; dat , the mistag 
probabilities, and the use of the effective discriminating variable. Since the physics 
parameters assumed in the three analyses are slightly different, the D~ results are adjusted 
to the more recent set of physics parameters listed in Table [I] before averaging. The 



combined amplitude plot is displayed in Fig. [10] and the corresponding numerical values 
are listed in Table 0. All values of Am s below 9.6 ps _1 are excluded at 95% CL. The 
combined sensitivity is 10.6 ps _1 . 

As the statistical correlation between this analysis and the previous ALEPH dilepton 
and lepton-kaon analyses |3], |5[] is very large, no significant improvement in sensitivity is 
expected if these latter analyses were included in the combination. 



12 Conclusion 

From a sample of 33023 inclusive lepton events, all values of Am s below 9.5 ps -1 are 
excluded at 95% CL using the amplitude method. This analysis supersedes the previous 
ALEPH inclusive lepton analysis |Q and provides the highest sensitivity and highest 
95% CL lower limit on Am s of any B° mixing analysis published to date |l|, 0, [| f|, |22|, |23| . 
Taking into account correlated systematic uncertainties the combination with the 
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ALEPH Combination 



A 95% CL limit 9.6 ps 
-0- sensitivity 10.6 ps -1 



data ± 1.645 a 

data ± 1.645 a (stat only) 




Ain (ps *) 



Figure 10: Measured Bg oscillation amplitude as a function of Am s for the combination of this 
analysis with the ALEPH D~ based analyses. 



ALEPH D~ based analyses yields Am s > 9.6 ps -1 at 95% CL. 
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Table 7: Combined measurements of the Bg oscillation amplitude A as a function of Am s (in 
ps _1 ), together with the statistical uncertainty cr^ at and the total systematic uncertainty (r s J st . 









Am s 


A±o%*±oJ( Bt 




Am s 




0.00 


+0.03±0.08±0.18 




7.00 


-0.15±0.30±0.18 




14.00 


±1.21±0.86±0.47 


1.00 


+0.16 ±0.11 ±0.18 




8.00 


-0.24 ±0.35 ±0.21 




15.00 


±1.98±0.99±0.54 


2.00 


-0.13±0.13±0.17 




9.00 


-0.05 ±0.40 ±0.23 




16.00 


±2.76±1.16±0.57 


3.00 


+0.13±0.16±0.19 




10.00 


±0.30 ±0.46 ±0.27 




17.00 


±2.86±1.37±0.61 


4.00 


+0.00±0.18±0.16 




11.00 


±0.37±0.54±0.37 




18.00 


±2.22 ±1.61 ±0.77 


5.00 


-0.10±0.22±0.18 




12.00 


±0.47±0.64±0.39 




19.00 


±1.85±1.88±0.98 


6.00 


-0.20±0.25±0.17 




13.00 


±0.65 ±0.75 ±0.42 




20.00 


±2.02±2.19±1.29 



References 

[1] ALEPH Collaboration, "Study of the B° — 5 S ° oscillation frequency using D~ — £ + 
combinations in Z decays", Phys. Lett. B 377 (1996) 205. 

[2] ALEPH Collaboration, "Study of B° s oscillations and lifetime using fully 
reconstructed D~ decays", Eur. Phys. J. C 4 (1998) 367. 

[3] ALEPH Collaboration, "An investigation of B° d and B° oscillation", Phys. Lett. B 
322 (1994) 441. 

[4] ALEPH Collaboration, "Limit on B° oscillation using a jet charge method", Phys. 
Lett. B 356 (1995) 409. 

[5] ALEPH Collaboration, "Combined limit on the B° s oscillation frequency", 
contribution PA08-020 to the 28 th International Conference on High Energy Physics, 
Warsaw, Poland, July 1996. 

[6] ALEPH Collaboration, "ALEPH: a detector for electron-positron annihilations at 
LEP", Nucl. Instrum. Methods A 294 (1990) 121, erratum 303 (1991) 393; 
B. Mours et al, Nucl. Instrum. Methods A 379 (1996) 101. 

[7] ALEPH Collaboration, "Performance of the ALEPH detector at LEP", Nucl. 
Instrum. Methods A 360 (1995) 481. 

[8] ALEPH Collaboration, "Improved measurements of electroweak parameters from Z 
decays into fermion pairs", Z. Phys. C 53 (1992) 1. 

[9] T. Sjostrand and M. Bengtsson, Computer Phys. Commun. 43 (1987) 367. 

[10] ALEPH Collaboration, "Heavy quark tagging with leptons in the ALEPH detector", 
Nucl. Instrum. Methods A 346 (1994) 461. 

[11] JADE Collaboration, "Experimental investigation of the energy dependence of the 
strong coupling strength", Phys. Lett. B 213 (1988) 235. 



22 



[12] ALEPH Collaboration, "A precise measurement ofT(Z — > bb)/T(Z — > hadrons)", 
Phys. Lett. B 313 (1993) 535. 

[13] O. Schneider, "Heavy quark spectroscopy, lifetimes and oscillations"", CERN- 
PPE/97-143, to appear in the proceedings of the 18 th International Symposium 
on Lepton-Photon Interactions, Hamburg, Germany, July 1997. 

[14] The LEP Collaborations ALEPH, DELPHI, L3, OPAL, the LEP Electroweak 
Working Group and the SLD Heavy Flavour Group, "A Combination of Preliminary 
Electroweak Measurements and Constraints on the Standard Model", CERN- 
PPE/97-154. 

[15] ALEPH, DELPHI, L3, OPAL Collaborations, "Combining heavy flavour electroweak 
measurements at LEP", Nucl. Instrum. Methods A 378 (1996) 101. 

[16] ALEPH Collaboration, "Measurement of mean lifetime and branching fractions of 
b -hadrons decaying to J/ if;", Phys. Lett. B 295 (1992) 396. 

[17] ALEPH Collaboration, "Measurement of the B° lifetime", Phys. Lett. B 322 (1994) 
275. 

[18] O. Leroy, "Etude du melange B° — B° dans V experience ALEPH au LEP", Ph.D. 
thesis, CPPM-T-1998-01, Universite de la Mediterranee, 1998. 

[19] H.-G. Moser and A. Roussarie, Nucl. Instrum. Methods A 384 (1997) 491. 

[20] ALEPH Collaboration, "Measurement of the b forward-backward asymmetry and 
mixing using high-p± leptons", Phys. lett. B 384 (1996) 414. 

[21] M. Beneke, G. Buchalla and I. Dunietz, Phys. Rev. D 54 (1996) 4419. 

[22] DELPHI Collaboration, "Search for B° s - B° oscillations" , Phys. Lett. B 414 (1997) 
382. 

[23] OPAL Collaboration, "An updated study of B meson oscillations using dilepton 
events", Z. Phys. C 76 (1997) 417; "A study of B meson oscillations using hadronic 
Z decays containing leptons", Z. Phys. C 76 (1997) 401. 



23 



